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1 INTRODUCTION 

1 —8 

None  of  the  recent  papers  from  the  RAE  , or  elsewhere  (eg  Refs  9 and  10), 
on  the  formulation  of  the  equations  of  motion  of  a deformable  aircraft,  can  be 
said  to  have  been  received  enthusiastically  or  even  given  general  satisfaction. 
There  has  been  criticism  for  incomprehensibility  and  for  naivety;  notation  and 
nomenclature  have  sometimes  been  off-putting;  and  there  have  been  particular 
criticisms  of  such  things  as  the  representation  of  the  perturbations,  the  choice 
of  frames  of  reference  etc.  What  is  suggested  here  is  a general  approach  which, 
it  is  believed,  is  simple  to  follow  and  can  be  particularized  to  suit  the  whim 
or  fancy  of  the  majority.  Thus  one  can  take  the  principal  frame  of  reference  to 
be  body-fixed  axes*,  mean-body  axes*,  constant-velocity  axes*,  or  a non-inertial 
frame  satisfying  some  other  (or  none  at  all)  conditions. 

We  will  restrict  ourselves  to  small  perturbations  of  a datum  motion  of 
constant  linear  velocity  and  zero  angular  velocity,  to  flight  in  a uniform  atmos- 
phere over  a 'flat  earth',  and  to  flight  of  an  aircraft  of  constant  mass  which  is 
such  that  no  part  of  it  moves  relative  to  any  other  part  during  the  datum  motion. 
Refinements,  such  as  allowing  for  the  effect  of  engine  rotation,  taking  other 
datum  motions  during  which  there  is  no  deformation,  etc,  can  be  made  on  the  lines 
developed  elsewhere  (Refs  6 and  8,  for  example). 

The  approach  is  briefly  as  follows.  We  take  as  a basic  frame  of  reference 
constant-velocity  axes*  which  are  such  that  they  coincide,  during  the  datum  motion, 
with  some  axes,  fixed  in  the  aircraft  during  the  datum  motion,  such  as  body-path 
axes  or  body-fixed  axes  (perhaps  the  principal  axes  of  inertia  ) . The  perturb- 
ation is  then  represented  by  a perturbation  with  no  deformation  which  moves  a set 
of  body-fixed  axes  from  coincidence  with  the  constant-velocity  axes  to  an  'arbit- 
rary' position  occupied  by  axes  which  we  call  no-deformation-body- fixed  axes, 
followed  by  deformations  which  are  such  that  the  translational  displacement  of 
any  particle  relative  to  the  no-deformation-body-fixed  axes  and  referred  to  those 
axes  is  precisely  a linear  function  of  a set  of  generalised  coordinates.  The 
equations  of  motion  are  then  obtained  using  Lagrange's  equation  for  a non-inertial 
frame,  to  obtain  the  equations  for  the  def ormational  freedoms;  and  applying  the 
principles  of  linear  and  angular  momentum  in  the  directions  of  and  about  the  no- 
deformation-body-fixed axes,  to  obtain  the  rigid  body  equations. 

Some  may  say  that  "we  want  to  treat  all  the  perturbations  in  the  same  way, 
we  don't  see  any  point  in  distinguishing  between  rigid  body  perturbations  and 

* These  are  defined  in  the  Glossary  at  the  end  of  the  paper. 
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deformations".  The  answer  is  that  you  do  not  have  to.  You  merely  take  the  particu- 
lar case  where  the  no-deformation-body-fixed  axes  coincide  with  the  constant- 
velocity  axes.  That  is,  put  the  specific  rigid  body  perturbations  to  zero  and 
include  in  the  deformational  freedoms  modes  which  are  exactly  (for  the  trans- 
lations) , or  a sufficiently  good  approximation  to  (for  the  rotations) , rigid 
body  perturbations. 

The  advantages  of  this  approach  are: 

(i)  The  principal  frame  of  reference  - the  no-deformation-body  fixed  axes 
can  be  made  any  non-inertial  frame*  that  one  desires,  or  left,  in  a sense, 
arbitrary,  by  a suitable  choice  of  deformation  modes.  Thus,  for  example,  modes 
satisfying  certain  conditions  will  make  it  body-fixed  axes,  while  modes  satisfying 
certain  other  conditions  will  make  it  mean-body  axes.  Alternatively,  as  already 
described,  we  can  make  the  principal  frame  of  reference  an  inertial  frame  - 
constant-velocity  axes. 

(ii)  The  generalised  applied  forces  that  appear  in  Lagrange's  equation  are 
obtained  by  calculating  the  virtual  work  done  in  infinitesimal  displacements 
assuming  the  frame  of  reference  (the  non-inertial  frame  which  is  our  principal 
frame  of  reference)  remains  stationary**.  With  our  suggested  approach  the 
contribution  from  any  applied  force  to  the  first  order  terms  (in  the  generalised 
coordinates)  in  the  generalised  forces  is,  since  the  deformation  is  precisely 
linear  in  the  generalised  coordinates,  given  merely  by  the  virtual  work  done  by 
the  first  order  perturbations  in  the  components  of  the  applied  force  referred  to 
the  no-deformation-body-fixed  axes.  Previous  developments  {eg  Ref  2)  have  been 
criticised  because  the  required  first  order  terms  in  the  generalised  forces 
involved  the  unperturbed  value  of  the  applied  force  - this  did  some  significant 
virtual  work  since  the  relevant  deformations  were  not  precisely  linear  in  the 
generalised  coordinates.  The  same  criticism  could  have  been  made  in  our  present 
development  if  we  had  applied  Lagrange's  equation  for  an  inertial  frame  (the 
constant-velocity  axes)  without  making  the  specific  rigid  body  rotations  all  zero. 

(iii)  The  overall  applied  forces  that  appear  in  the  other  equations  of  motion  - 
the  rigid  body  equations  - are  also  such  that  there  is  no  contribution,  from  the 
unperturbed  value  of  an  applied  force,  to  their  linear  constituents  (coefficients 
of  the  generalised  coordinates) . 

* As  already  stated  it  can  be  made  an  inertial  frame  by  omission  of  the  specific 
body  freedoms. 

**  That  is  that  the  position  of  the  frame  of  reference  is  not  affected  by  the 
infinitesimal  displacements. 
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In  this  paper  the  development  of  the  equations  of  motion  is  made  entirely 
in  a dimensional  form.  There  is  however  little  difficulty  in  writing  them  in  a 
non-dimensional  form,  as,  for  example,  was  done  in  Ref  3;  though  some  would 
question  whether  anything  is  gained  by  so  doing  rather  than  working  with  the 
dimensional  form  and  a system  of  units*  appropriate  to  the  problem  being 
considered. 

2 THE  SEMI-RIGID  MODEL 

Let  the  constant-velocity  frame  of  reference  be  identified  by  constant- 

velocity  axes  0 x y z . The  location  of  the  no-deformation-body-fixed  axes 
c c c c 

0 x y z can  then  be  specified  by  the  translations  and  rotations  which,  when 
n n'n  n 

applied  to  the  constant-velocity  axes,  would  make  the  two  sets  of  axes  coincident. 
Let  these  translations  and  rotations  be,  in  succession,: 

(c)  (c)  (c) 

(i)  Translations  Xj  , yj  , Zj  in  the  directions  of  the  respective 
constant-velocity  axes. 

(ii)  Successive  rotations  i|»,  6,  <)>  about  the  carried  axes  0,0,0  . 

z y x 

These  translations  and  rotations  will  be  taken  as  the  six  specific  'body-freedom' 
generalised  coordinates**.  They  uniquely  specify  the  position  of  the  no- 
deformation-body-fixed axes.  Following  these  rigid  body  perturbations  we  assume 
the  aircraft  is  deformed  so  that  the  position  of  a particle,  relative  to  the 
origin  of  the  no-deformation-body-fixed  axes,  is  given  precisely  by 
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are  funct: 

ates  (x^,  y^,  z^)  of  the  unperturbed  particle,  and  S is  the  axes  transformation 
matrix  (see  equation  (2-4)  below). 

Writing  the  linear  velocity  components  of  the  constant-velocity  axes, 
referred  to  those  axes,  as  (u^,  v^,  w^.},  and  remembering  that  their  angular 


* Such  as  the  aero-normalised  system  proposed  by  Hopkin 
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**  However,  if  desired,  as  stated  in  the  Introduction,  these  steps  just  described 
can  be  omitted  and  good  approximations  to  the  body-freedoms  included  amongst 
the  def ormational  freedoms. 
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velocity  is  zero,  it  is  easily  seen  (Ref  1,  Appendix  A)  that  the  linear  and  angular 
velocities  of  the  no-deformation-body-fixed  axes  are,  referred  to  those  axes. 
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where*  S is  the  axes  transformation  matrix 


and 


S = I - A + . . . 
<P 


(2-2) 


(2-3) 


(2-4) 

(2-5) 


in  which,  for  any  column  vector  {<p  6 ip} 


and 


0 - ip  6 

\p  0 - <p: 

0 <p  0! 

0 0 0 

0 0 -ip 

0 (p  0 


(2-6) 


(2-7) 


Consequently  the  linear  velocity  of  a particle,  referred  to  the  no-deformation 
body-fixed  axes,  is  given  by  (of  Ref  6,  for  example). 


* See  Ref  14,  Appendix  M and  Ref  1,  Appendix  A;  and  note  that  S is  orthogonal 
ST  - s-'. 
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where  D is  Che  operator  d/d t , and,  in  addition  to  (2-3)  to  (2-5),  we  have  made 
use  of  the  fact  that 
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(2-8a) 


for  any  vectors  { d>  9 i|t},  {x  y z}. 


2 . 1 Various  particularisations 

It  has  been  pointed  out  in  Ref  6 that  if  the  modal  matrix  R satisfies  the 
conditions 


y~^3mAx  R = 0 


(2-9) 


where  6m  is  the  mass  of  an  element  and  the  summation  is  for  all  the  elements  of 
the  aircraft,  then  the  no-def ormation-body-f ixed  axes  are  identical  with  the  mean- 
body  axes  (satisfying  the  condition  that  the  kinetic  energy  relative  to  the  axes 
is  a minimum).  This  condition  (2-9)  means  that  the  deformation  modes  are  ortho- 
gonal to  the  specific  body  freedoms  with  respect  to  the  mass  distribution. 

If  in  addition  to  (2-9)  we  impose  the  condition 


1 1 


Zsm 
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(2-10) 
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which  makes  the  origin  of  the  no-deformation-body-fixed  axes  at  the  aircraft  CG 
when  there  is  no  deformation,  then  from  (2-1),  the  origin  of  these  axes  (alias  the 
mean-body  axes)  stays  at  the  CG  whatever  the  deformation.  Equation  (2-10)  means 
that  the  specific  body  freedoms  are  mutually  orthogonal  with  respect  to  the  mass 
distribution. 


6 . 


xf  = yf  = zf 


Mentioned  in  the  same  paper  is  the  case  when  at  the  particle 
0 


R = 0 

has  null  second  and  third  rows 


SR 

Sx, 


SR 

sy* 


has  null  third  and  first  rows 


SR 

5zr 


has  null  first  and  second  rows. 

J 


(2-11) 


The  no-deformation-body-fixed  axes  are  then  identical  with  body-fixed  axes  fixed 
in  position  and  orientation,  in  the  aircraft,  at  their  origin. 


Another  case  of  interest  occurs  when  we  make  the  specific  body  freedom 
(c) 


coordinates  Xj''  , 4>  etc)  all  zero  and  put 


R = 


I -A 


(2-12) 


The  last  three  of  these  modes  are  then  good  approximations  to  the  rigid  body 
rotations,  and  the  previous  three  are  rigid  body  translations.  Moreover  the 
no-deformation-body-fixed  axes  are  then  identical  with  the  constant-velocity  axes. 
The  approach  then  becomes  that  favoured  by  many  flutter  analysts  (of  Ref  7). 

With  this  representation  the  coordinates  of  a particle  relative  to  the  comstant- 
velocity  axes  is 
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which  compares  with  our  basic  representation 


(c)  T 

X.  + % Xf  + R 


xf  I + R q i + xi 


- A + 

xf 


(2-14) 


Thus  the  two  representations  agree  for  small  perturbations  if,  in  (2-13),  we 
replace  n by  n + 6,  and  then  put 


(2-15) 


On  occasion  it  may  be  desirable  or  convenient  to  specify  the  deformation  of 
the  aircraft  in  terms  of  translations  and  rotations  of  strips  of  the  aircraft  and 
of  rotations  of  'flap'  parts  of  the  strips  relative  to  the  main  parts  of  their 
own  strips.  It  will  be  seen  by  comparison  with  the  analysis  of  Ref  5,  or  other- 
wise, that  a good  approximation  to  such  a representation  is  achieved  by  taking 


T 1 

S.  A Qa  F.  + S 

*if  Xsif  *if  1 « 


(2-16) 


at  the  ith  strip. 
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In  this  expression  the  <t>.^  etc  express  the  orientation  of  some  strip-attached 

axes  relative  to  the  no-deformation-body-fixed  axes  in  the  datum  state*;  the 

x ^ etc  are  the  datum  coordinates  of  a particle  in  the  s ;rip-attached  frame  of 

reference;  x . , = x,  . (y  . c = z . . = 0)  is  the  hinge  position;  and  K. , F. , f. 

sif  nsi  •'sif  sif  ’ i*  1 1 

are  respectively  translational  rotational  (approximate)  and  flap  rotational 
(approximate)  modal  matrices.  For  our  present  purpose  we  will  simplify  the  above 
expression  by  assuming  that  in  the  datum  state  all  the  strip-attached  axes  have 
the  same  orientation  as  the  no-deformation-body-fixed  axes  in  the  datum  state. 
Thus  equation  (2-16)  becomes  then 


R 


iXsif  ^Ssi 


(2-17) 


at  the  ith  strip.  This  is  an  extreme  simplification.  A more  representative  model 
would  be  to  take 
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(2-18) 


where  the  second  option  would  be  used  for  strips  such  as  those  on  a fin.  With  a 
modal  matrix  of  the  form  (2-17),  since  in  this  case 
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* The  strip-attached  axes  are  attached  to  the  strip  at  the  strip  reference  point, 
which  is  in  the  main  part,  not  the  flap  part,  of  the  strip.  During  the  datum 
motion  they  are  fixed  in  the  strip.  During  perturbed  motion,  according  to  the 
above  modal  matrix  R (equation  (2-16)),  the  deformation  of  the  strip  will 
only  differ,  from  a deformation  which  is  purely  a rotation  of  the  flap  part  of 
the  strip,  by  terms  which  are  of  second  order  in  the  generalised  coordinates. 

Thus,  when  the  perturbations  are  small,  we  can  say  that  the  strip-attached  axes 
are  almost  fixed  in  the  main  part  of  the  strip.  We  are  assuming,  by  definition.  Oil 

that  the  flap  angle  is  zero  in  the  datum  state.  Ref  5 is  more  general  being 
written  in  terms  of  datum  values  of  the  flap  angle. 


iMMl 


I 


where  (x^,  y^,  zif^  ’-s  t*le  str^P  reference  point  in  the  datum  state,  we  can 
include  good  approximations  to  the  rigid  body  modes,  as  in  (2-12),  by  taking 


K.  = K.  I -A 

t 1 > 


(2-20) 


'i  ■ [?i  0 D 


(2-21) 


q . [?j  o o] 


(2-22) 


In  the  more  general  representation  (2-16),  since  then 


x-  f + S x . . 

if  6 . . Slf 

if  1 


(2-23) 


one  would  have,  instead  of  (2-21),  to  take 


?.  = F. 

l l 


0 q*1 

♦if  ♦if 


THE  APPLIED  FORCES 


(2-24) 


Let  a typical  local  applied  force,  referred  to  the  no-deformation-body-fixed 
axes,  be  (the  bar  suprascript  denotes  typical) 


= e-  + e.  . . . e q,  + e e e x. 
f 1 n 1 x y z 1 


+ e e e.  <(>  + ... 

<f>  9 1 p 


(3-1) 


where  the  coefficients  e j , e^,  e^  etc  may  be  differential  operators.  In  the 
following  two  sub-sections  we  give  the  contributions,  from  a set  of  such  forces. 
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! 


to  the  overall  and  generalised  applied  forces.  Then  in  sections  4 to  6 we  obtain 
these  contributions  when  applied  forces  are,  respectively,  aerodynamic  gravit- 
ational and  structural.  Since  this  paper  is  intentionally  short,  propulsive 
forces  have  not  been  considered  separately  - they  were  considered  in  Refs  I to  6 - 
it  being  left  to  the  reader  to  do  that,  or,  alternatively,  to  include  them  with 
the  aerodynamic  forces. 


For  the  particular  case  (equation  (2-12))  mentioned  in  section  2.1,  when 
n is  increased  to  n + $ and  R put  equal  to  R , equation  (3-1)  becomes 


to  the  fact  that  the  no-deformation-body-fixed  axes  are  differently  oriented,  when 
the  first  of  these  is  zero  and  not  the  second,  compared  with  the  case  when  the 
second  is  zero  and  not  the  first*.  Of  course  when  one  takes  the  above  form  one 
will  normally  make  all  the  specific  body  freedom  coordinates  (x|C\  9 etc)  zero. 

* In  this  latter  case  the  no-deformation-body-fixed  axes  coincide  with  the 

constant-velocity  axes,  while  in  general  they  do  not.  The  coefficients  of 

q , to  q , in  (3-2)  were  therefore  found  from  the  condition 
^n+1  ^n+6 


<P 


9 

4> 


(3-2a) 


Oil 
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3 . 1 The  generalised  applied  forces 


By  considering  the  virtual  work  done  by  a set  of  forces  (3-1)  in  infinite- 
simal displacements,  assuming  the  no-deformation-body-fixed  axes  to  be  stationary* 
(of  equation  (2-1)),  their  contributions  to  the  generalised  applied  forces  are 
easily  seen  to  be 
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(3-3) 


3.2  The  overall  applied  forces 

The  contributions  from  the  typical  set  of  applied  forces  (3-1)  to  the  trans- 
lational force  on  the  aircraft,  and  moment  about  the  origin  of  the  no-deformation- 
body-fixed  axes,  referred  to  those  axes,  are  respectively 


~x(n)* 

■ y 

- (n) 

e 

y(n) 

i— i 

f (n) 

z(n) 

-(n) 

g 

(3-4) 


That  is  that  the  position  of  the  frame  of  reference  is  not  affected  by  the 
infinitesimal  displacements. 
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and 


n 


M(n) 

n 


n 


(3-5) 


In  the  particular  case,  (equation  (2-12)),  mentioned  in  section  2.1,  when  n is 
increased  to  n + 6 and  R put  equal  to  R , the  above  two  expressions  become, 
using  (3-2)  instead  of  (3-1) 


(3-6) 


0 


+ . . . 


(3-7) 


4 THE  AERODYNAMIC  FORCES 

The  boundary  condition  to  be  satisfied  by  the  air  motion  at  the  surface  of 
the  body  involves  the  velocity  of  the  surface  normal  to  itself.  The  velocity 
vector  of  a particle  is  given  by  equation  (2-8)  in  terms  of  its  components 
referred  to  the  no-deformation-body-fixed  axes.  The  slope  vector  of  the  body 
surface  relative  to  the  same  axes  will  be  merely  a function  of  the  generalised 
coordinates  qj  to  q^  (of  equation  (2-1)).  The  local  aerodynamic  force  vector 
will  therefore  be  a function  of  the  three  quantities 
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The  coefficients  X3| , etc  in  the  above  two  equations,  (4-4)  and  (4-5),  may 
be  differential  operators. 

For  the  particular  case  (equation  (2-12))  mentioned  in  section  2.1,  when  n 
is  increased  to  n + 6,  R is  put  equal  to  R , and  the  specific  body  freedom 
coordinates  made  zero,  then  the  expression  for  the  local  aerodynamic  force  vector 
has  the  form  (of  equations  (4-2)  and  (3-2)) 
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e , + e . . . . e . q,  + 

af  al  a,n+6  1 


where 


’af  *al 


ea,n+l  ea,n+3 


(4-6) 


a,n+l 

a,n+l 


e»  e*  e«  D 
x y z 


(4-7) 


0 0 
a,n+4  . . . a,n+6 


a,n+4 

®a,n+4 


e . e . e . A 
x y z u^ 


er  e*  er  D - A 

4>  8 \p  e 


(4-8) 


The  generalised  forces  are  then  given  by  the  obvious  modification  of  (4-3)  - 
(c) 

x|  , <J>  etc  zero,  the  use  of  (4-6),  and  the  replacement  of  R by  [ R I - AXf ] . 
4. 1 The  use  of  strip  aerodynamics 


If  the  modes  have  been  expressed  in  a strip  form,  such  as  (2-17),  then  the 
user  may  find  it  more  convenient  or  more  understandable  to  have  the  aerodynamic 
forces  expressed  in  terms  of  the  overall  force,  moment  etc  on  a strip.  Assuming 
two-dimensional  aerodynamics  for  each  strip,  and  no  inter-strip  interference*, 
we  can  write'’  the  overall  forces  on  a strip,  when  in  the  unperturbed  state  the 
strip-attached  axes  have  the  same  orientation  as  the  constant-velocity  axes,  as 


X.c  + X..  X..  X.i  X..  ur'  - u,  + 

if  IX  1Z  10  10  1 f 

0 0000  wfs)  - wc 

l f 

Z.r  z..  z..  z..  Z.t  qfs) 

if  lx  iz  10  i6  ^i 


(4-9) 


where  the  superscript  (s)  indcates  the  resolution  along  the  strip-attached  axes; 
ujS\  wjs)  are  components  of  the  velocity  of  the  strip  reference  point;  and 


* Not  even  between  two  strips  where  one  lies  aft  of  the  other.  For  the  purpose 
of  this  paper  things  are  being  kept  simple. 
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The  deformation  of  the  strip,  other  than  the  flap  rotation,  will  only  affect 
higher  order  terms  in  (4-9),  since,  as  we  have  pointed  out,  the  assumed  behaviour 
(equations  (2-1)  and  (2-17))  only  differs  from  that  of  a rigid  aerofoil  and  a 
rigid  flap  by  a term  which  is  of  second  order  of  smallness.  Thus,  from  (4-9), 
using  (4-11),  the  overall  forces  on  a strip  resolved  along  the  no-deformation- 
body-fixed  axes  are 
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Similarly  the  overall  moment  on  the  strip,  about  the  strip  reference  point,  and 
the  hinge  moment,  referred  to  the  no-deformation-body-fixed  axes  have  the  forms 
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where  (of  equations  (2-1)  and  (2-17). 
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(4-21) 


and  the  " superscript  indicates  a summation  over  the  flap  part  of  the  strip. 


Carrying  out  the  comparisons  enables  us  to  express  various  stripwise  summations. 


of  the  coefficients  that  appear  in  the  expression  (equation  (4-2))  for  the  local 
aerodynamic  force  vector,  in  terms  of  the  strip  aerodynamic  coefficients  etc. 

These  summations  appear  in  the  expressions  for  the  overall  and  generalised  aero- 
dynamic forces  (equations  (4-3),  (4-4)  and  (4-5)),  bearing  in  mind  that  they  relate 


to  the  case  when  the  modal  matrix  R has  the  form  (2-17).  We  will  not  detail  the 


relationships  here,  nor  will  we,  in  the  present  paper,  give  any  further  consider- 
ation to  the  use  of  strip  theory  aerodynamics. 


5 THE  GRAVITATIONAL  FORCES 

On  each  particle  of  the  aircraft  there  will  be,  according  to  the  standard 
assumption,  a gravitational  force  6mg  acting  in  the  direction  (vertically 

downwards)  of  the  normal  earth-fixed  axes*.  If  the  orientation  transformation 
matrix  which  takes  us  from  these  axes  to  the  constant-velocity  axes  is  , 


where 


(5-1) 


is  similar  to  S ; then  the  matrix  which  performs  the  orientation  transformation 
from  the  normal  earth-fixed  axes  to  the  no-deformation-body-fixed  axes  is  SS* 

f 

Thus  the  local  gravitational  force  vector  referred  to  the  no-deformation-body- 
fixed  axes  is 


* See  Glossary  of  nomenclature. 
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where  is  the  third  column  of 

f f 

Consequently,  from  the  formulae  of  section  3,  the  gravitational  contributions  to 
the  generalised  and  overall  applied  forces  are,  respectively. 


and 


(5-3) 


(5-4) 


(5-5) 


For  the  particular  case,  mentioned  in  section  2.1  (equation  2-12)),  when  n 
is  increased  to  n + 6,  R put  equal  to  R , and  the  specific  body  freedom 
coordinates  made  zero,  then  the  expression  for  the  local  gravitational  force 
vector  has  the  form  (of  equations  (5-1)  and  (3-2)) 


23 


e 

f 

g 


(nr 

g 

(n) 

g 

(n) 

g 


6mgJlcj) 


The  generalised  forces  then  are 


6 THE  STRUCTURAL  FORCES 


(5-6) 


(5-7) 


Resulting  from  such  things  as  the  stresses  in  the  material,  friction  between 
different  components  of  the  structure  etc,  there  will  be  a force  on  any  particle 
of  the  aircraft  which  can  be  called  a structural  force.  The  meaning  of  this  term 
is  considered  more  fully  in  Appendix  A of  Ref  2.  This  structural  force  vector 

will  depend  only  on  the  deformation  of  the  aircraft  and  so  we  can  write  its 

/ (n)  (n)  (n 

components  [e  , f , g 
\ s s s 

the  form 


referred  to  the  no-deformation-body-fixed  axes,  in 
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where  the  coefficients,  eg^  etc,  may  be  differential  operators.  This  approxi- 
mates to  any  contribution  not  resulting  from  a stress  by  a contribution  of  the 
above  form. 

The  structure  cannot  exert  any  overall  force  or  moment  on  itself  and  so  it 
is  easily  seen  (of  Appendix  A of  Ref  2)  that  the  following  conditions  must  be 
satisfied : 


01  1 


Consequently  from  the  formulae  of  section  3 there  is  no  structural  contribution 
to  the  overall  force  and  moment,  while  the  contribution  to  the  generalised  forces 
is 
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For  the  particular  case,  mentioned  in  section  2.1  (equation  (2-12)),  when 
n is  increased  to  n + 6,  R put  equal  to  R , and  the  specific  body  freedom 
coordinates  made  zero,  then  the  expression  for  the  local  stiffness  force  vector 
becomes  ( of  equations  (6-1)  and  (3-2)) 
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where 
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The  generalised  forces  (structural  contribution)  then  are 
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where  use  has  been  made  of  equations  (6-2)  to  (6-5) . Many  readers  will  be  amazed 
to  see  non-zero  terms  outside  the  (11)  submatrix  in  this  equation  and  in  particular 
to  see  the  term  {^Ax^Aeg^.}  in  the  bottom  right  hand  corner.  "How  can  a structure 
have  an  elastic  stiffness  in  degrees  of  freedom  which  to  all  intents  and  purposes 
are  rigid  body  rotations?"  one  asks.  Consider  the  case  when  the  only  freedoms 
are  6,  . The  local  structural  force  vector  is  then 
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referred  to  the  no-deformation-body-fixed  axes  and 
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referred  to  the  constant-velocity  axes.  The  coordinates  of  a typical  particle 
referred  to  the  cons tant-ve loci ty  axes  are  similarly 
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where  B,Q  is  the  second  order  term  in  the  expansion  of  S (of  equation  (2-4)). 
q>u 

Thus  the  moment  of  the  force  distribution  (6-12)  about  the  origin  of  the  constant- 
velocity  axes  is 
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and  both  these  terms  are  zero  as  a consequency  of  (6-3)  . Also  the  virtual  work 
done  by  (6-12)  in  infinitesimal  displacements  6$,  69,  Sip  is,  using  a result  from 
Ref  2, equation  (57), 
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(7-3) 


(7-4) 


(7-5) 


Q.  is  the  generalised  applied  force  obtained  by  the  method  of 
virtual  work  - in  the  assumed  virtual  displacement  the  frame 
of  reference  is  regarded  as  stationary 


and  q.  is  the  generalised  coordinate  of  the  ith  degree  of  freedom 
where  the  freedom  is  one  such  that  the  position  and  orient- 
ation of  the  frame  of  reference  is  independent  of  it. 

The  expression  on  the  left  hand  side  of  the  equation  (7-1)  is  known  as  the 
effective  force  for  the  degree  of  freedom  q^  . Using  the  expressions  for  the 
velocities  of  the  no-deformation-body-fixed  axes  (equations  (2-2)  and  (2-3))  and  Oil 

for  the  location  of  a particle  relative  to  those  axes  (equation  (2-1))  we  then  find 
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that  the  column  vector  of  the  reversed  effective  forces  for  the  deformational 
freedoms  is 
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Applying  the  principles  of  linear  and  angular  momentum  respectively  in  the 
directions  of  and  about  stationary  axes  which  are  instantaneously  coincident  with 
the  no-deformation-body-fixed  axes  we  obtain  the  relationships 
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where  X 


(n)  v(n)  (n) 


are  the  resolutes  along  the  no-deformation-body-fixed 
(n)  w(n)  „(n) 


axes  of  the  total  applied  force;  and  L^“'  , M^“' , are  t^le  comPotlents  about 

those  axes  of  the  total  applied  torque  about  their  origin.  Thus  the  effective 
forces  for  the  body  freedom  equations  are  given  by  the  right  hand  sides  of  the 
above  equations  and  are  (using  equations  (2-1),  (2-3)  and  (2.8)) 

(i)  for  the  translational  equations 
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(ii)  for  the  rotational  freedoms 
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8 THE  EQUATIONS  OF  MOTION 

The  equations  of  motion,  which  can  be  described  as  the  equation  of  the 
applied  and  effective  forces,  have  already  been  stated  in  general  terms  in 
section  7 (equations  (7-1),  (7-7)  and  (7-8)).  These  must  be  satisfied  whatever 
the  perturbation;  and  in  particular  when  the  perturbation  of  the  datum  motion  is 
zero,  and  when  that  perturbation  is  small. 


8. 1 The  equation  of  equilibrium 


When  the  perturbations  are  zero  the  effective  forces  are  zero  and  so  adding 
together  the  various  contributions  to  the  generalised  and  overall  forces  we  find 
that  we  must  have 
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8.2  The  perturbation  motion  equations 

When  we  collect  together,  from  sections  4 to  7 the  various  contributions 
to  the  applied  and  effective,  generalised  and  overall  forces,  we  find  that  for 


small  perturbations  the  following  equation  must  be  satisfied,  in  addition  to 

(8-0, 


where  the  various  constituent  matrices  are: 
(i)  Ponderous  inertia 


(ii)  Gravitational 
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(iii)  Structural 
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(iv)  Aerodynamics 
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8.3  Various  particularisations  of  the  perturbation  motion  equations 

Of  the  cases  mentioned  in  section  2.1,  probably  the  one  of  most  interest  is 
that  when  the  modal  matrix  has  the  form  (2-12),  with  n increased  to  n + 6 and 
R put  equal  to  R , and  the  specific  body  freedom  coordinates  are  made  zero.  The 
perturbation  motion  equation  then  is 


«>ii^  * [«yl  * - feijH 


1n+6 


(8-9) 


where  the  ponderous  inertia  matrix  is  given  by  (8-5),  and  the  other  matrices 


01  1 


are 
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(ii)  Gravitational  (of  equation  (5-7)) 

[G..]  = 0 


(iii)  Structural  (of  equation  (6-10)) 


(8-12) 


When  one  compares  this  form  of  the  equations  of  motion  with  that  of  the  previous 
section  one  sees  that  there  are  differences  in  the  individual  contributions  to 
the  (13),  (23),  (31)  and  (33)  submatrices.  However  as  a consequence  of  the 
equilibrium  equation  (equation  (8-2))  the  total  difference  in  the  second  of  these 
submatrices  is  clearly  zero.  In  the  (33)  matrix  the  total  difference  is 


L 
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A + A 
e e 

, af  sf 


(8-13) 


However  in  the  datum  state  the  total  force  on  any  particle  must  be  zero*  and  so 


SragS-rj, 


' M ' [‘I' ] 


at  every  particle 


(8-14) 


and  so  the  above  quantity  (equation  (8-13))  is  also  zero.  Equation  (8-14)  also 


makes 


{Mv  * IrT(\£  * \f)  ■ 0 


(8-15) 


and  so  the  differences  in  the  (13)  and  (31)  submatrices  are  also  zero. 

Finally,  in  connection  with  this  particularisation,  some  may  think  that  the 
appearance  of  datum  values  of  the  applied  forces  in  the  various  expressions  for 
the  coefficients  (equations  (8-11)  and  (8-12))  is  a denial  of  one  of  the  advantages 
claimed  (af  section  1)  for  the  present  approach.  However  this  is  just  a conse- 
quence of  writing  the  typical  coefficient  matrix 


en+4  en+5  en+6 


in  terms  of  the  coefficients  appearing  in  the  general  form  (af  equations  (3-2) 
and  (3-1)). 

The  above  is  the  form  the  equations  of  motion  take  when  the  principal 

frame  of  reference  used  is  a set  of  constant-velocity  axes.  There  have  however 
. 11-13 

been  quite  a number  of  people  who  have  advocated  using  mean-body  axes  as  the 

principal  frame.  This  means  merely  taking  a set  of  deformation  modes  which 


* The  equilibrium  equations  of  section  (8.1)  are  only  a truly  sufficient  state- 
ment if  the  semi-rigid  model  is  exact.  If  (8-14)  is  satisfied  then  (8-1)  to 
(8-3)  are  automatically  satisfied  but  the  converse  is  not  necessarily  true. 


satisfy  (2-9).  The  consequent  effect  on  the  general  perturbation  equations  of 
section  8.2  is  that  the  inertial  and  gravitational  coupling  between  the  deform- 
ational  and  rigid  body  freedoms  is  eliminated.  With  the  particular  mean-body  axes 
which  have  their  origin  at  the  aircraft  CG  then  (of  equation  (2-10),  the  inertial 
cross-coupling  between  the  translational  and  rotational  body  freedoms,  and  the 
direct  and  cross  gravitational  terms  for  the  rotational  body  freedoms,  are  also 
eliminated. 

In  section  2. 1 the  expression  of  the  modal  behaviour  in  terms  of  the  dis- 
placements etc  of  strips,  and  the  way  the  modes  could  be  chosen  to  make  the  no- 
deformation-body-fixed  axes  always  body-fixed,  were  also  discussed.  These 
particularisations  however  make  no  change  to  the  perturbation  motion  equation  of 
section  8.2  and  it  is  merely  a matter  of  making  the  appropriate  substitution  for 
the  modal  matrix  R . 
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37 


At  times  the  ability  of  the  structural  dynamicist's  approach  to  predict  low 
frequency,  predominantly  rigid,  perturbations  from  the  aircraft's  datum  flight 
state  has  been  questioned.  This  doubt  is  clearly  a consequence  of  the  fact  that 
the  flutter  investigator,  when  considering  symmetric  perturbations,  normally  does 
not  include  any  freedom  in  fore  and  aft  translation,  and  so  the  lowest  frequency 
characteristic  solution  that  he  obtains  is  the  short-period.  We  will  show  that 
the  much  lower  frequency  phugoid  can  also  be  obtained,  by  considering  its  deter- 
mination for  the  rigid  aircraft,  using  the  same  approach. 

A. 1 The  structural  dynamicist  takes  constant-velocity  axes  as  his  principal 
frame  of  reference  and  so  the  perturbation  motion  equation  for  the  'rigid'  air- 
craft is  that  given  at  the  beginning  of  section  8.3  with  n = 0.  Thus  the  various 
matrices  are: 


(i)  Gravitational 


0 


(ii)  Structural 
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0 


I 


A A 
Xf  esf 


(A-  1 ) 


(A- 2) 


(iii)  Aerodynamic 
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(iv)  Inertia 


As  mentioned  in  section  3,  propulsive  forces  can  be  treated  separately  rather  than 
include  them  in  the  aerodynamic  forces.  For  our  present  purpose  it  is  convenient 
so  to  do.  We  will  assume  that  the  local  propulsive  force  has  constant  components, 
referred  to  axes  fixed  in  the  aircraft,  during  perfectly  rigid  perturbations  from 
the  datum  motion.  This  means  that  in  the  expression  for  the  local  propulsive 
force  (of  the  expression  for  the  typical  local  force  equation  (3-2))  the  coef- 
ficient matrices 


— 

e e e 

and 

e , e „ e . 

px  py  pz 

p<f>  p0  p*p 

f 

f 

px 

p <t> 

2 

°px 

°p4>  J 

are  both  zero.  Consequently,  by  comparison  with  the  derivation  of  the  aerodynamic 
matrix,  or  otherwise,  we  find  that  we  have  a propulsive  matrix 


0 


A 

f SpfJ 


(A-5) 


w 


and  the  perturbation  motion  equation  is 
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{[Aij]»2  * [G.j]  * [E..]  * [Pjj]  - [Qij]}  f-fC> 


(A-6) 


where  we  have  now  replaced  q^  to  q^  by  symbols  which  indicate  the  nature  of 
the  perturbations*. 


A. 2 For  the  determination  of  the  phugoid  we  need  only  consider  symmetric  perturb- 


ations of  unbanked  (*f>  = 0)  flight,  and  so  putting  y(c  , <£*,  \p*  all  equal  to 


zero  in  (A-6)  it  becomes 


mET  - X-D  - X* D ( > 6mzjD2  - X:D  + X.w 
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6 x f z f 
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0 0 X , + X , 

af  pf 


o o - ^|(gaf  ♦ gsf  + gpf)Zf  ♦ (eaf  ♦ es£  + epf)x£f 


= 0 (A- 7) 


Now  from  the  equations  of  equilibrium  ((8-2)  and  (8-3)  with  the  propulsive  forces 


included)  we  have,  since  <f'£  = 0, 


if  + Xpf  = m8  Sin  ®f 


Z , + Z , 
af  pf 


mg  cos  O 


(A-8) 


In  addition,  as  pointed  out  in  section  8.3,  the  total  force  on  any  particle  must 
be  zero,  and  so  (af  equation  (8-14)) 


* It  is  hoped  that  this  will  cause  no  confusion  since  the  specific  body  freedom 
coordinates  were  deleted  at  a somewhat  earlier  stage  (of  section  8.3).  The 
asterisk  has  been  attached  to  the  rotational  symbols  to  indicate  that  these  do 
not  precisely  represent  rotations. 
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e,  + e,  + e,  = 6mg  sin  ft, 
af  sf  pf  6 f 


gaf  + 8sf  + 8Pf 


- 6mg  cos 


(A-9) 


Therefore 


(e  , + e , + e , )x, 
af  sf  pf  f 


<8af  + 8sf  + 8pf>Zf 


■ (Es“‘)e 

■ • (E-*) 


g sin  ft. 


g cos  ft 


(A- 10) 


A. 3 It  is  known  from  experience  that  during  the  phugoid  the  incidence  of  the 

aircraft  remains  approximately  constant  throughout.  We  will  therefore  assume 

that  the  perturbations  are  such  that  the  incidence  remains  constant.  With  the 

perturbation  9*  any  finite  straight  line  on  the  aircraft,  in  a y^11^  = 

constant  plane,  remains  a straight  line  but  it  is  rotated  through  an  angle 

tan  !0*  about  the  axis  as  well  as  undergoing  some  translation  and  slight 

change  in  length.  The  incidence  is  defined  by  the  direction,  of  the  velocity  of 

some  reference  point  on  the  aircraft,  relative  to  such  a line.  We  will  take  as 

the  reference  point  the  particle  (x,  = y,  = zf  = 0) . Then  since,  from  equation 

t t r t 

(2-8),  for  our  present  case  , for  this  particle 


we  require 


u(n)  - w(n)9* 
m ra 

W(n)  + u(n)e* 
m ra 


, .(c)  1 

uf  + X1  I 


Wf  + 21 


(wf + z5c))< 


uf  + x - + 2 7 

wf + 2ic)  + (uf + xic)v* 


(A- 11) 


(A-  1 2) 


+ One  has  to  go  through  the  steps 

Jc)  _ , 


0 etc 
n + 6 
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te 


• (c)  • (c) 

wfX]  - UfZj 

2 2 -(c)  .(c) 

uf  wf  + ufXj  + wfz] 


A .(c) 

,2  X1 


2 i<c)  + 
.2  Z1 


(A- 13) 


For  level  flight,  in  which  case 


= Vf  cos  0^ 


V sin  0^ 


(A- 14) 


this  condition  becomes 


•(c) 

0*  = sin  0 


.(c) 


f "vT  cos  0f  + 


(A- 15) 


Putting  this  condition  into  the  first  two  equations  of  (A-7),  neglecting  the 
third  equation,  assuming  that  in  the  datum  state  the  aircraft  CG  coincides  with 
the  particle  x^  **  yf  = zf  = 0,  and  making  use  of  (A-8)  and  (A— 14),  gives 


(A- 16) 
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B 


For  our  present  purpose  it  will  be  good  enough  to  assume  that  the  two  6 coef- 
ficients, Xf  and  Z*  , are  zero  and  that  the  other  aerodynamic  coefficients  are 
o 0 

constant.  We  then  see  that  the  two  non-trivial  solutions  of  (A-16)  have  a 
circular  frequency  w given  by 


+ X*  sin  0,  cos  ©„  + 
x f £ 


2£ 

vf 


s 0,  + Z • sin  0, 


0-28 


vf  si"2°fJ 


- |(X.  cos^Q,  - Z.  sin^©. 


The  temporal  behaviour  of  the  solution  is  cos  (cut  + constant) 

where  p = ~ jx*  cos2©,  + (X*  + Z*)  sin  ©,  cos  0,  + Z*  sin2©,}  . (A- 18) 

2m  ' x t zx  f fz  f 


To  quantify  the  above  solution  let  us  take  a simple  aerodynamic  model  which 

accords  with  the  assumptions  that  we  have  already  made.  Thus  we  assume  that  the 

aerodynamic  force  consists  of  (i)  a force  acting  normal  to  the  direction  of  motion 

of  the  reference  particle  (xf  = yf  = zf  = 0)  which  is  a function  of  the  velocity 

of  that  particle  and  the  incidence,  and  (ii)  a force  acting  in  the  direction  of 

motion  of  the  reference  particle  which  is  a function  of  the  incidence  times  the 
2 

(velocity)  of  that  particle.  In  particular  we  write  the  former  as 


l\  3Nf  3Nf  I 

N = V0jNf  + 3V7  (V0  - V + (a0  ' V[ 


(A— 19) 


and  the  latter  as  -D 


where 
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and  (of  equations  (A— 1 1 ) and  (A— 12)) 
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The  components  of  the  overall  aerodynamic  force  in  the  directions  of  the  constant- 
velocity  axes  then  are 

• (c)  • (c) 
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(ao  ' V 


+Uf3Vj(V0-V+  •••  * 


and  so,  since  (of  equations  (3-6)  and  (4-4)) 


Appendix 


45 


With  these  values  it  can  be  confirmed  that  the  coefficients  of  6 in  (A-25)  and 
(A-29)  also  agree,  as,  of  course,  they  should. 

A. 5 We  can  proceed  to  evaluate  the  frequency  and  decay  rate  of  the  phugoid  from 
equations  (A- 17)  and  (A- 18)  when  u^ , w^  are  given  by  (A- 14),  we  find  that 


X.  sin  0 + X.  sin  ©,  cos  © = - 2V. 

z f x f f f 


sin  ©j  cos  ©^  + sin2©^ 


v2  8Nf  • 2© 
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(A-36) 
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x f z f f 
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(A- 37) 
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X*  cos  0,  - Z.  sin  ©.  + (X.  - Z.)  sin  0,  cos  0_ 
X fz  f zx  f f 


= ~ 2V^|d^ ^cos2©^  - sin2©^j  + N^(2  sin 


+ N^(2  sin  ©^  cos  0^)^ 


2 3Nf 

- Vf  (2  sin  ©f  cos  ©f) 


(A-38) 


2 2 2 3Nf 

X.  sin  ©,  + (X.  - Z-)  sin  0,  cos  0,  - Z.  cos  ©f  = - 2V£Nr  - V7  rrA  (A-39) 

z fxz  f fx  r fff  9V^ 


2 2 
X.  cos  ©f  + (X.  + Z.)  sin  ©f  cos  ©f  + Z.  sin  ©f  = - 2VCD 

X I Z X L ^ Z I f f 


(A-40) 


Consequently,  from  (A- 17)  and  (A- 18) 
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Bearing  in  the  mind  that  we  have  (from  (A-8) , (A--14),  (A-30)  and  (A— 3 1 )) 


X , sin  6.  - Z r cos  - mg 
pf f pf f 


(A-43) 


it  is  more  revealing  to  rewrite  the  expression  for  the  frequency  in  the  form 


2 - 


!f  3Nf 

mg 


- ("*  * 5)  ' (£) 


9V, 


(A-44) 


This  reduces  to  Lanchester's  classical  solution  (w  = (g/V^)/2  ; y = 0)  when  the 
propulsive  force  is  zero  and  is  zero;  for,  similar  to  (A-43). 

X , cos  0.  + Z c sin  0. 

D = ^ l (A-45) 

Vf 


and  so  this,  and  y , are  zero. 

If  only  3N^^3V^  and  (N^  + (mg/v|))  are  zero  the  above  solution  ((A-44)  and 
(A-42))  is  the  'better  approximation'  obtained  by  Etkin*  (Ref  12,  p 331).  Without 
either  of  these  assumptions  we  have  an  even  better  approximation  to  the  true 
solution;  and  it  is  clear  that  as  accurate  a determination  of  the  phugoid  can  be 
obtained  by  the  present  approach  as  has  been  obtained  by  any  other.  We  made  a 
number  of  simplifying  assumptions  - constant  incidence,  no  genuine  deformational 
freedoms,  simple  aerodynamic  model  and  so  on  - as  our  purpose  was  merely  to 
demonstrate  the  adequacy  of  the  suggested  approach  rather  than  make  an  accurate 
calculation  for  some  particular  case.  When  one  is  concerned  with  this  latter 
aspect  the  above  assumptions  can  be  rejected  as  desired  and  indeed  in  numerical 
work  there  is  probably  little  to  be  gained  by  making  most  of  them. 


* He  writes  the  damping  ratio  and  frequency  in  terms  of  the  derivative  Cfy  of 

the  thrust  coefficient  with  respect  to  the  (airspeed/unperturbed  airspeed),  but 
this  is  a little  deceptive  for  he  has  already  assumed  (pages  earlier)  that  the 
thrust  is  constant  and  so  C-pv  is  simply  related  to  the  thrust  coefficient  C^,  : 

CTV  " “ 2CT 
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GLOSSARY  OF  NOMENCLATURE 


(i)  Frames  of  reference 
Body-fixed  axes 

Body-path  axes 

Constant-velocity  axes 

Mean-body  axes 

No-deformation-body-f ixed 
axes 


Normal  earth-fixed  axes 

Principal  axes  of  inertia 

Strip-attached  axes 

(ii)  Miscellaneous 

Aircraft  reference  point 

Axes  transformation  matrix 

Strip 

Strip  reference  point 


(all  right-handed  orthogonal  cartesian) 

axes  whose  origin  and  orientation  are  fixed 
in  a small  material  portion  of  the  body 

axes  whose  origin  is  at  a material  point  of 
the  body  and  such  that  the  velocity  of  that 
point  is  in  the  direction  of  the  x axis. 

axes  having  constant  linear  and  angular 
velocity  relative  to  an  inertial  frame.  (The 
ones  used  have  zero  angular  velocity.) 

axes  orientated  so  that  the  kinetic  energy 
relative  to  the  axes  is  a minimum. 

an  arbitrary  non-inertial  frame.  Any  perturb- 
ation can  be  visualised  as  a rigid  body 
perturbation  which  moves  the  body-fixed  axes 
from  coincidence  with  the  constant-velocity 
axes  to  coincidence  with  these  axes,  followed 
by  a deformation. 

axes  fixed  relative  to  the  earth  with  the 
z-axis  vertically  downwards. 

axes  with  origin  at  the  centre  of  gravity  and 
such  that  the  three  products  of  inertia  about 
the  axes  are  zero. 

axes  with  origin  at  a material  point  of  a 
strip  (the  Strip  reference  point)  and  such 
that,  in  our  strip  model,  the  deformation  of 
the  main  part  of  the  strip  relative  to  these 
axes  is  of  second  order  of  smallness. 


the  material  point  on  the  aircraft  which  is 
the  origin  of  the  body-fixed  axes. 

the  matrix  which  by  pre-multiplying  a column 
vector  changes  the  axes  directions  in  respect 
to  which  the  vector  is  resolved  into 
components . 

a narrow  slice  of  the  aircraft  lying  more  or 
less  for  and  aft  and  consisting  of  one  or 
two  almost  rigid  portions. 

a chosen  material  point  on  the  main  part  of 
a strip  taken  as  the  origin  of  the  strip- 
attached  axes. 
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LIST  OF  SYMBOLS 

skew-symmetric  matrices  involving  <p , 0,  ip  etc  (see 

equation  (2-6)) 

ponderous  inertia  coefficient 

aerodynamic  hinge  moment  on  ith  strip 

aerodynamic  hinge  moment  coefficients  for  ith  strip 
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-E. 

l 

E.  . 
ij 

F. 

l 

-G. 

l 

G. 

l 

G.  . 
ij 

I 

J. 

l 


M 


M. 


, M*  etc 

x $ 

N 


N. 

l 


differential  operator  d/dt 
generalised  structural  force 
structural  stiffness  coefficient 

rotational  (approximate)  modal  matrix  for  ith  strip 
generalised  gravitational  force 
gyrostatic  force  (of  equation  (7-4)) 
gravitational  stiffness  coefficient 
unit  matrix  (3  x 3) 

a certain  coupling  force  between  the  rotational  body  freedoms 
and  the  def ormational  freedoms  in  Lagrange's  equations 
referred  to  a non-inertial  frame  (see  equation  (7-5)) 
a certain  matrix  formed  from  the  elements  of  {$  0 i|»}  - see 
equation  (2-7) 

translational  model  matrix  for  the  ith  strip 

rolling  moment 
rolling  moment  on  ith  strip 

< 

rolling  moment  coefficient  (aerodynamic) 
aerodynamic  rolling  moment  coefficients 
pitching  moment 

(pitching  moment  on  ith  strip 
pitching  moment  coefficient  (aerodynamic) 
aerodynamic  pitching  moment  coefficients 
yawing  moment 

(yawing  moment  on  ith  strip 
yawing  moment  coefficient  (aerodynamic) 
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LIST  OF  SYMBOLS  (continued) 


Q; 


Qij 


s,  s, 


4> 


W 

X,  Y,  Z 


x. , Y.  , 

z. 

1 1 

l 

x. . X* 

etc 

x 4> 

Y..  Y* 

etc 

x $ 

Z. . Z‘ 

etc 

x 

e,  f,  g 

e . , e , 

e . 

1 X 

4> 

f.,  f , 

fA> 

1 X 

4> 

g.  » g » 

V 

1 X 

> etc 


e* , e* 

x $ 

f.,  f; 
x’  <t> 

8x’  gi 


1 

8 

m 

6m 

n 

P.  q.  r 
9; 


etc 


u,  v,  w 
x,  y,  z 


generalised  force 

matrix  relating  angular  velocities  and  orientation  = 

I 0 - sin  9 

0 cos  <j>  sin  <J>  cos  0 

0 - sin  $ cos  $ cos  9 

aerodynamic  coefficient 
modal  matrix  (see  equation  (2-1)) 

axes  transformation  matrix.  Absence  of  a subscript  means  the 
arguments  are  <(> , 9 » 'i' . S=I-A^  + B^q  ... 
centrifugal  potential  function 

kinetic  energy  relative  to  the  frame  of  reference 
overall  force  resolutes 

(overall  force  resolutes  on  ith  strip 
overall  force  resolute  coefficients  (aerodynamic) 

aerodynamic  force  resolute  coefficients  - they  may  be 
differential  operators 

components  of  local  force  vector 

local  force  vector  coefficients  - they  may  be  differential 
operators 


local  aerodynamic  force  vector  coefficients  - they  may  be 
differential  operators 

flap  modal  vector  at  ith  strip 
acceleration  due  to  gravity 
third  column  of 
mass  of  the  aircraft 
mass  of  a particle 

number  of  deformational  degrees  of  freedom 
angular  velocity  resolutes 
generalised  coordinate 
time 

linear  velocity  resolutes 
particle  position  resolutes 
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I*  "1 


4> 


'I'. 


>»  0,  *P 


■if’  9if’  *if 


Dressings 


LIST  OF  SYMBOLS  (continued) 


position  resolutes  for  origin  of  no-deformation-body-fixed  axes 
relative  to  the  origin  of  the  constant-velocity  axes 


angle  of  inclination 
in  datum  motion 

angle  of  bank  in  datum  ^ 
motion 

nose-azimuth  angle  in 
datum  motion 


orientation  angles  of  constant- 
velocity  axes  relative  to  normal 
earth-fixed  axes 


angle  of  flap  rotation  (approximate) 

orientation  angles  of  no-deformation-body-fixed  axes  relative 
to  constant-velocity  axes 

orientation  angles  of  the  strip-attached  axes,  relative  to  the 
no-deforraation-body-f ixed  axes,  in  the  datum  state 


(i)  Subscripts 

A quantity  is  relative  to  or  about  the  origin  of  a particular  set  of  axes 

where 


c denotes  the  constant-velocity  axes 

h denotes  axes  with  the  same  orientation  as  the  strip-attached  axes  but  with 
origin  at  the  strip  hinge 
n denotes  the  no-deformation-body-fixed  axes 

s denotes  the  strip-attached  axes 

a indicates  aerodynamic 

1 indicates  the  value  during  the  datum  motion 
g indicates  gravitational 

i indicates  the  ith  strip 

s indicates  structural 

(ii)  Superscripts 

Bracketed  superscripts  denote  the  axes  of  resolution  with  the  same  signific- 
ance as  in  (i)  above,  also 

" indicates  a summation  over  the  flap  part  of  a strip 
T denotes  the  transpose  of  a matrix 
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LIST  OF  SYMBOLS  (concluded) 


• (dot)  denotes  derivative  with  respect  to  time 

~ (tilde)  is  of  ephemeral  significance 

(bar)  denotes  typical 
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